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We extend the Polyakov-Nambu-Jona-Lasinio (PNJL) model for two degenerate flavours by in- 
cluding the effect of the SU(3) measure with a Vandermonde (VdM) term. This ensures that the 
Polyakov loop always remains in the domain [0,1]. The pressure, energy density, specific heat, speed 
of sound and conformal measure show small or negligible effects from this term. However various 
. . . quark number and isospin susceptibilities are all found to approach their respective ideal gas limits 

00 ■ around 2 T^. We compare our methods with other similar approaches in PNJL model and also 

, present a quantitative comparison with Lattice QCD data. 
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00 ; I. INTRODUCTION 

I— I. Recently, there is a lot of interest in the studies of thermodynamics of strongly interacting matter using the 
i-C ' Polyakov loop enhanced Nambu-Jona-Lasinio (PNJL) model [1-3]. This model couples the chiral and deconfinement 
^p-i, order parameters through a simple-minded coupling of the NJL model [4] with the Polyakov loop model [5] . The two 
^ • major thrusts in recent times have been to estimate various thermodynamic observables using this model (see e.g. 
D [ [6-10]), and to make systematic improvements of the model [11-13]. Another set of important result has come from 
j-^^' similar studies in chiral quark models that go beyond the mean field treatment [14]. 

In this note we deal with the improvement of the Polyakov loop model and describe some of its consequences, 
. remaining within the domain of mean field analysis. The Polyakov loop model used in many of the recent literature 
I> ' is the one given in Ref.[3]. The Polyakov loop <& has been treated here as a Z(3) spin field [15]. Using this model we 
, estimated [6] a very sensitive observable - the quark number susceptibility (QNS) and also the higher order coefficients 
in the Taylor expansion of pressure in quark number chemical potential /xq. Comparison with the data from Lattice 
' QCD (LQCD) [16] showed that the QNS in the PNJL model and LQCD agree quite well both qualitatively and 
quantitatively. The fourth order coefficient C4 showed qualitative agreement but had a quantitative difference at high 
temperatures. Some of us further extended the PNJL model to include isospin chemical potential ni [8]. The isospin 
number susceptibility (INS) and its derivative with respect to no and fj,i were obtained. In this case the fourth order 
derivative C4 was quite consistent with lattice data, but the INS was not. A possible reason for such departures is 
that the mean-field treatment of the PNJL model is insufficient. But then it should have affected the coefficients 
systematically i.e., all the fourth order coefficients should deviate further from LQCD data than the second order 
coefficients. 

There are however other simpler reasons that should be considered first. The PNJL model is only a model which 
d ] can mimic some of the characteristics of a fundamental theory like QCD and its discretized version LQCD. Moreover, 
the parameters like the couplings and masses are quite different in the PNJL model and the LQCD simulations. Thus 
some quantitative difference is naturally expected. Apart from these we made an important observation in [8] that 
$ has a big role to play in the behaviour of these coefficients. We pointed out how the quantitative differences could 
be caused by the behaviour of $ as a function of temperature and chemical potentials. The most important physical 
problem in the simple-minded PNJL model is the following. $ being the normalized trace of the Wilson line L, which 
is an SU(3) matrix, should lie in the range < 4> < 1. But it was found to be greater than 1 at temperatures above 
2Tc (see Fig. 2 in Ref. [8]). The natural way to cure this problem is to consider a proper Jacobian of transformation 
from the matrix valued field L to the complex valued field $ which will then constrain the value of $ to ^ < 1. This 
is quite a well known construction in SU(N) matrix model (see e.g. [17-19]), in certain variations of Polyakov loop 
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model ([11, 20]), as well as in QCD motivated phcnomenological models (see [21] and references therein). Also this is 
ubiquitous in various strong coupling effective theories of Lattice QCD (see e.g. [23]). 

Here we introduce the Vandermonde term in the Polyakov loop model in a conceptually different way than that in 
the earlier models. In the next section we discuss our approach. In section III we show the changes in measurements 
of the susceptibilities and various other quantities due to the VdM term. The final section contains our conclusions. 



II. FORMALISM 

At a temperature T, the SU(3) Wilson line is given by L(x) = 'Pexp{ig J^^^ AQ{x.)XadT), where g is the gauge 
coupling, Aq (a = 1,2, ...8) are the time-like components of the gluon field, Aa are the Gell-Mann matrices and r 
the imaginary time in the Euclidian field theory. The Polyakov loop is defined as 4> = trL/3 and its conjugate is 
i" = trL^/3. Since L is itself a SU(3) matrix so $ < 1. The gluon thermodynamics can be described as an effective 
theory of the Polyakov loops [5] . On the other hand quark thermodynamics can be effectively described in terms of 
NJL model [4], and the two are coupled to obtain the PNJL model (e.g., [3]). The thermodynamic potential in this 
model can be obtained in terms of the sigma and pion condensates and the thermal average of the Polyakov loop. 

However the version of the PNJL model [3] leads to <!> > 1 for T > 2Tc. To rectify this anomaly, the authors of 
Ref. [3] have recently proposed a complete modification of the Polyakov loop model [11], motivated from the strong 
coupling results used by Fukushima [2] Our aim in this work is also similar, but the approach is somewhat different. 
We retain the Polyakov loop potential of [3, 6, 8] but treat it as a matrix model. Also the way we define pressure is 
quite different as discussed below. 

We first outline our scheme using an arbitrary matrix model for the Wilson line L, which for simplicity is assumed 
to be a potential V[L] = V[$,i']. In the following equation, we express the partition function for this theory first as 
a path integral over L and then over the fields <f> and 

PLe-*^[*^*l = /'n^L(x)e-^^[*'*] (la) 

= J[$(x),#(x)]d$(x)di(x)e-+^[*'*1 (lb) 

where, 2?L is the SU(3) Haar measure, J[<I>, <t] is the Jacobian of transformation (also called Vandermonde determinant, 
see e.g. Ref. [24]) from L to ($, $), and is given as J[$, l>] ee (27/247r2)(l - 6 l>$ + 4 ($3 + $3) _ 3 ($$)2-)^ Qur 
interest then would be to obtain the pressure which is given by, 

d InZ /dV\ 1 , , 

where, v denotes the physical volume of the system and () denotes thermal averaging. The last approximation holds 
in the infinite volume limit. 

The role of the Jacobian is to be understood as follows. First, it is a factor reweighting the field configurations and 
hence significantly affects all thermal averages. However the Jacobian is not explicitly space-time dependent, there is 
no extra term to be averaged in Eqn. 2 as one might expect when redefining the path integration from L to $ (Eqn. 
la to Eqn. lb). A typical example of such a dependence would be if we were considering say a Fourier transform 
of the fields. In case of a free field this kind of dependence of the Jacobian on the volume and temperature is very 
important in obtaining the correct partition function. 

Thus, in our mean field treatment we have to carefully incorporate the effect of the Jacobian and this is the main 
aim of this paper. The effect of the Jacobian is refiected in the mean fields ($) and ($), and we express the pressure 
as, 

p^--vm,{^)). (3) 

V 

To relate to pure glue theory, we now replace the potential density V/v by a Landau-Ginzburg type functional U, 
given by [3], 

U(<!>,^,T) b2(T)- ba , ^ bi .2 

W ^ =-^^^~-|(^' + $^) + ^(a>$)' , (4) 
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with 



&2 (T) = ao + ai ( ^ 



a2 \ — j +0-3 1 Y 



(5) 





FIG. 1: $ and P/Psb for k = 0(ro = 0.27 GeV) and k = Q.5{To = 0.2555 GeV). The value of is 0.270 GeV. 

To make a saddle point approximation to the mean fields, the potential density U was minimized w.r.t. ^ and $ in 
Rcf. [3]. These were then used to obtain pressure P = —U. The coefficients (i=0,l,2,3) and hj (j=2,3,4) were 
fitted from Lattice data of pressure in pure gauge theory, and Tg is precisely the transition temperature Tc = 270 
MeV [25-27]. As T oo, PjT^ 167r^/90. However, to take care of the effect of the Jacobian as discussed above, 
we now propose to minimize the following modified potential, 



Z^($,$) 



y4 



y4 



Kln[J($,$)] 



(6) 



where k is a dimensionless parameter to be determined phenomenologically. The mean field value of pressure is still 
obtained from the relation P — —U. A very simple example of this approach is demonstrated in the appendix. Note 
that the Jacobian term is considered as an extra effective term in the modified potential density implying a sort of 
normalized volume factor. This is quite natural as the form of Eqn. lb implies that there is a Jacobian sitting at 
each and every space-time coordinate, depending on the value of <f> and $. 




T/T, 




FIG. 2: $ and PjPsB for k = 0(To = 0.27 GeV) and k = {0.22T^ /T^){To = 0.2555 GeV). Here Tc = 0.270 GeV. 



With the new minimization condition all the coefficients should be estimated afresh. Instead, we retain the values 
of tti and bj obtained in [3] and tune only the values of Tq and k. This is equivalent to a correlated modification of 
the Qi and bj keeping Tq fixed at 270 MeV. 
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We show the variation of the Polyakov loop and the pressure P normalized to Stefan-Boltzmann (SB) pressure 
PsB for pure gauge theory, as a function of temperature. In Fig.l we have used a small non-zero constant value of 
K = 0.05. In Fig. 2 we find similar behaviour for a temperature dependent k = Q.22Tq /T^ . In both the figures the 
K = curves are for the Polyakov loop model without the VdM term. Thus the parameter space of k is quite open 
at this stage. 

Within the range of temperatures (T < 3Tc) where the Polyakov loop model is supposed to be a good description 
of the system, our approach and that of Ref. [11] give similar results. The reason behind this is that one can suitably 
adjust the parameters in both approaches. However, our method for introducing the VdM potential as discussed 
above, is very much different from that of Ref. [11]. The main difference is that the pressure computed in Rcf. [11] 
includes the VdM term. Thus the coefficient of the VdM term requires an inverse temperature dependence, so that on 
a naive extrapolation to high temperatures, the pressure does not blow up with the logarithm of the Jacobian. In that 
case another problem crops up with the remaining part of the thermodynamic potential, which at high temperatures 
has no bound, contrary to the claim that $ ^ 1 as T — oo. Precisely because 3> should go to 1 as T ^ oo we believe 
that the VdM term should be very important at high temperatures to constrain the maximum value of <& to 1. 

The exercise for introducing a VdM term for the Polyakov loop model itself has nothing new to offer. Even without 
it the potential U was able to describe the pure glue theory quite well. However its importance becomes evident 
in the PNJL model. The Polyakov loop has a coupling to the fermionic part as will be seen in the corresponding 
thermodynamic potential below, which forces the $ to be greater than 1, and more so as the chemical potential is 
increased. The VdM term can inhibit such a behaviour. 

The thermodynamic potential of the PNJL model [3, 6, 8] is given as. 



n 



^ 2T / ^{ln[l + 3(a> + $c-(^/-''/)/^)e-(^/-^/)/^ + c-3(^/-''/)/^ 



+ ln[l + 3($ + $e-(^^+^/)/^)c-(^/+w)/7' + c-3(^/+«)/^]}- ^ 6 / -^Ef9 {A^ ~ p 



i2nr 



(7) 



Here quark condensates for the two light flavors u and d are given by (t„ ~< uu > and Od =< dd > respectively, and 
the respective chemical potentials are /i„ and ^d- Note that /io = {^u + fJ'd)/'2 and fij = {piu — fJ-d)/'^- The quasi-particle 

energies are Eu^d ~ \/P^ + '^u d' where m„_d = mo — ^Giau,d — ^G2<Jd,u are the constituent quark masses and toq 



is the current quark mass (we assume flavour degeneracy). Gi and G2 are the effective coupling strengths of a local, 
chiral symmetric four-point interaction. We take Gi = G2 = G/4, where G is the coupling used in Rcf. [3]. A is 
the 3-momentum cutoff in the NJL model. U ($, $,T) is the effective potential for $ and $ as given in Eqn. 4. We 
locate the transition temperature in this model from the peaks in the temperature variation of d^/dT and da^.d/dT. 
Similar to the case of the Polyakov loop model we would now obtain the mean fields by minimizing. 



^ = ^-Kln[J($,<I)] 



(8) 



The coefficient k in the VdM term can in general have some temperature and/or chemical potential dependence. 
Here we take a constant value k = 0.2 which suffices for the purpose of the present work. To set this value we looked 
at the two important quantities affected by the VdM term. First one is $ which decreases with the increase of k 
and hence decreases the pressure. Second one is the transition temperature which increases with k. Thus we try to 
optimize k to get both the pressure and the transition temperature as close as possible to the LQCD results for two 
quark flavours. 

On a naive extrapolation of this model to large chemical potentials, the $ and $ should grow towards 1 (decon- 
fincment at large chemical potential) even at very low temperatures. Thus again the logarithmic term blows up. So 
if pressure is computed including the VdM term as is done in Ref. [11], an anomalous logarithmic divergence would 
come up. There may be some new physics that can obscure such terms by making k as /i 00. But that would 
again run into a problem in restricting $ in the domain < < 1 • 

Apart from the difference in the treatment of the VdM term we would now remove the condition $ = $ used in 
[11], since it has important implications for susceptibilities. 

Before going over to our results let us take a digression to the Lattice computation of <I>. On the lattice $ is 
computed from the relation [15], 



$(r) = exp(- AF,,-(oo,r)/2T), 



(9) 
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where, Afg^(oo, T) = Fqqipo, T) — Foq{T), and Fqq{r, T) is the free energy of a pair of heavy quark and anti-quark at 
a separation r at a temperature T. This has been used to define a renormahzed Polyakov loop in lattice simulations 
of both pure gluon [28, 29] and full QCD [30]. In fact the data of [28] was used to obtain the different parameters 
of the Polyakov loop model in [3], and is being used by us here, and in that sense <& is the renormalizcd Polyakov 
loop. But even in this exercise the <I> in the Polyakov loop model of [3] goes to 1 at large T and is thus different 
from lattice results for T > Tc- On the lattice the value of <I> goes above 1 for T > Tc- It has been argued that 
since the $ measured in lattice simulations is a renormalizcd quantity, it is no more a character of the group SU(3) 
and is thus not limited to values below 1. From Eqn.9, it is evident that <& > 1 only when AFgg(oo,r) < 0, and 
this can be very easily seen to be true in the lattice simulations and happens for T > Tc- Now, the free energy 
Fqq(r,T) can be considered to be composed of three components, namely, a confining potential, a screening potential 
and an entropy part. For low temperatures the confining part is dominant and AFqq{oo,T) > 0. In the deconfincd 
phase for large distances, the screening potential drops out so the entropy part is dominant which could lead to 
AFqq{oo,T) ~ -T A Sqq{T) < 0, whcrc ASqq{T) = Sqq{T) - Sqo, and Sqq{T) denotes the entropy of the system 
with a pair of quark and anti-quark. However the heavy quarks as such are not expected to contribute significantly 
to the entropy and it seems natural to have ASqq{T) = 0, and thus AFqq{oo,T) = for T > Tc. Instead the value 
is negative on the lattice and AFqq{oo,T) — > — oo as T — > oo, leading to $ — > oo. One has to then worry about 
what can bend it down towards 1 at asymptotic temperatures as was observed by Cava and Jcngo in perturbative 
evaluation of $ [31]. However this perturbative calculation also points to the fact that as the temperature is lowered 
from asymptotic values the $ is greater than 1. Also recent continuum estimates in chiral quark models [32] using 
dimensional reduction find close agreement with both lattice and perturbative calculations. 

On the other hand another lattice computation of the Polyakov loop in pure glue theory uses a renormalization 
dependent on temperature instead on the lattice spacing and finds the values to remain below 1 at least upto T ~ 3.5Tc 
[33]. We thus admit that the state of affairs with the lattice computation of $ is not very clear to us at this stage. 
There is a missing link from quantum computations to our matrix model mean-field computations. 



III. RESULTS AND DISCUSSIONS 



A. PNJL Model: Pressure, specific heat and speed of sound 



Now wc discuss the results for the PNJL model with VdM term. Here the fl' as given in Eqn. 8 is minimized with 
respect to the fields and all the thermodynamic quantities are obtained using these values. The peaks of the d^/dT 
and d(Tu,d/dT curves, as shown in Fig. 3(a), differs by 5 MeV. Their average position, which is at 230 MeV, is taken 
as the transition (or crossover) temperature Tc- In spite of the significant difference of Tc in the PNJL model with 
the corresponding LQCD value of 192 (7) (4) MeV [34], the thermodynamic quantities when plotted against the scaled 
temperature T/Tc show similar behaviour. We shall henceforth show the temperature dependences in terms of T /Tc. 

As mentioned earlier we are using an optimized value of k = 0.2. The temperature dependence of the fields are 
shown in Fig. 3(b). It agrees reasonably with that of the LQCD results as shown in Fig. 1 of Ref. [35]. The scaled 
pressure P/Psb is plotted in Fig. 3(c). It shghtly overestimates the LQCD pressure [36]. However it agrees well with 
the recent LQCD results for 2-1-1 flavors with almost physical quark masses [37]. 

Now, the energy density e is obtained from the relation. 



d /Q 



dT \T 



= -T 



df 



n 



(10) 



The rate of change of energy density e with temperature at constant volume is the specific heat Cy which is given as, 



_ de 



= - T 



9T2 



(11) 



The square of velocity of sound at constant entropy S is given by. 
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(12) 



The conformal measure is given by. 
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FIG. 3: (a): Peaks in d<l>/dr and dcr/dT sets the Tc at around 230 MeV. (b): $ and cr as functions of T/Tc. 
Note: In this figure o — G{au + Od)- 



C = A/e 



A = g - 3P 



(13) 





FIG. 4: (a): Temperature dependence of energy density e and specific heat Cv- (b): Temperature dependence of squared speed 
of sound and conformal measure A/e. The arrows on the right show the corresponding SB hmit. 



These quantities are plotted in Fig. 4. At higher temperatures the Cy is slightly lower than the values obtained in 
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[6]. However, the velocity of sound and the conformal measure remain unaltered in the whole range of temperatures. 
Thus the VdM term affects Cy but not quantities involving ratios of pressure and energy density e.g. v1 and C. It 
is interesting to note that our earlier [6] as well as the present work, have been able to predict the value of v1 quite 
well when compared to the recent LQCD results [37]. We hope similar encouraging results would be obtained on the 
lattice for the specific heat. 



B. Taylor expansion of Pressure 
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FIG. 5: The Taylor expansion coefficients of pressure in quark number and isospin chemical potentials as functions of T/Tc. 
Symbols are LQCD data [16]. Arrows on the right indicate the corresponding ideal gas values. 

The Taylor expansion coefficients of pressure with respect to chemical potentials have been the focus of comparison 
of PNJL and LQCD results [6, 8, 11, 22]. Here we have expanded the scaled pressure (P/T^) in a Taylor scries for 
the quark number and isospin number chemical potentials, (Iq and /x/ respectively. 



2^4 



-EE 



n=0 j=0 



; k = n- j, 



where, 



1 9" {P{T,fio,fii)/T^) 



(14) 



(15) 



The n = odd terms vanish due to CP symmetry. Even for the n = even terms, due to flavour degeneracy all the 
coefficients c^^ with j and k both odd vanish identically. We evaluate all the 10 nonzero coefficients (including the 
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pressure at = = 0) upto order n = 6 and compare them to LQCD data. These coefficients were evaluated in 
[6, 8] and certain differences were found w.r.t. LQCD data. We shall now discuss the effects of the VdM term on 
these coefficients. 

The coefficients we deal with are given by, 



CniT) = - 



1 9"(P(T,^o)/T4 



Mo=0 



(16) 



ifrr., 1 a"(P(r,Mo,M/)/r4 



n> 1. 



(17) 



/iO=0,/J7=0 



We present the QNS, INS and their higher order derivatives with respect to /io in Fig. 5. We have plotted the LQCD 
data from Ref. [16] for quantitative comparison. At the second order (Fig. 5(a)) we find that the QNS C2 compares 
well with the LQCD data upto about 1.2 Tc- Thereafter the PNJL values rise up towards the SB limit, while the 
LQCD values saturate at about 80% of this limit. The INS C2 also shows similar behaviour, but at lower temperatures 
it goes slightly above the corresponding LQCD values. There is no significant difference of C2 with and without the 
VdM term. However C2 was close to the LQCD result without VdM term [8], but now at high temperatures it goes 
above the LQCD values and approaches C2. Thus at high temperatures these coefficients overestimate the LQCD 
results but both are almost equal to each other, similar to that observed on the Lattice. This was not so without the 
VdM term [8]. 

Now we discuss the 4*'' order coefficients (Fig. 5(b)). The values of C4 in the PNJL model with VdM term matches 
closely with those of LQCD data for the full range of temperatures. This is in contrast to that found without the 
VdM term [6] where they were close only upto T ^ l.lTc. The VdM term does not affect the coefficient C4 which 
agrees well with LQCD data for the full range of T. Also both these coefficients approach each other as well as the 
corresponding SB limit. At the 6*'' order (Fig. 5(c)) the coefficients do not seem to be affected by the VdM term. 

Thus we write down the salient features regarding the Taylor coefficients in this modified PNJL model: 

• All the coefficients start approaching their respective SB limit around 2Tc. 

• Both the QNS and INS approach each other at 2Tc. This is also true for their corresponding responses to quark 
chemical potential given by the 4*'' and 6*'' order coefficients. 

• At high temperatures, except C2 and Cj, all the coefficients compare well quantitatively with the LQCD data. 

• The main effect of the VdM term is to move C2 and C4 close to their respective SB limits. 
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FIG. 6; (a): $ (solid lines) decreases and $ (dotted lines) increases as a function of j-io/T {/ij = 0) at low temperatures and 
almost equal and constant at high temperatures, (b): $ (solid lines) and ^ (dotted lines) are equal and almost constant as a 
function of j-ii /T (^0 = 0). 
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We have emphasized the role of the Polyakov loop in obtaining the values of the Taylor cocfRcients in our earlier 
works [6, 8]. In those works we found firstly that the Polyakov loop goes above 1 at high temperatures and also has 
a significant dependence on /iq but not on jij. Here as shown in Fig. 6, the VdM term restricts the value of $ within 
1, and also the dependence at higher temperatures is almost negligible. Thus even the splitting between $ and $ 
has almost disappeared. We note here that though we let <& and $ to be different, they come out to be almost equal 
at high temperatures. This is in contrast to imposing $ = $ for the full range of temperatures as done in Ref. [11]. 
The difference between $ and $ is responsible for the difference of C2 and in the intermediate temperatures. 

To complete the comparison with the LQCD data we have looked at the flavour diagonal (cJJ") and flavour off- 
diagonal {dlf) susceptibilities defined as, 



uu _ _n ~^ 



iO ^("-2)2 



and 



4 ' "4 
The 2-nd order flavour diagonal and off-diagonal susceptibilities arc given by, 
Xuu{T, = 0, = 0) d'^P{T, fJ-d) 



T2 df,l 
Xud{T, flu = 0, = 0) d'^P{T, fJ-d) 



rp2 



P-lL=tJ-d=0 



P.u=tJ.d=0 



2c"''. 



(18) 



and 




These are shown in Fig. 7. Except Cj", all the other LQCD diagonal and off-diagonal coefficients are close to their 
respective ideal gas values from 1.2rc onwards. The most striking discrepancy without the VdM term w.r.t the LQCD 
data was (see [8]) in the 2-nd order flavour off-diagonal susceptibility d^'''- c!^"^ signifies the mixing of u and d quarks 
through the contribution of the two disconnected u and d quark loops. While the LQCD data shows that this kind 
of correlation between the u-d flavours are almost zero just away form Tc, the PNJL model results remained non-zero 
even upto 2Tc. Adding the VdM term this part of the PNJL physics is now consistent with LQCD results. Below 
1.2Tc there is still a large quantitative difference between the PNJL and LQCD results for c!^"^. Obviously the VdM 
term is not expected to affect the results at low temperatures significantly. At the moment it is not clear what physics 
lie behind the difference between PNJL and LQCD results for at high temperatures and 03'' at low temperatures. 
Perhaps the quark masses may hold an answer. 
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IV. SUMMARY 



In this work the PNJL model of Ref. [3, 6, 8] has been extended by introducmg a VdM term. The hiiportant change it 
brings about is to set the upper hmit of the Polyakov loop to 1. With this model we have studied some thermodynamic 
properties of strongly interacting matter with the light flavours u and d within a certain range of temperature T, and 
small values of chemical potentials fiQ and /ij. In principle the VdM term affects all thermodynamic quantities. We 
adjusted the parameters in the model so that the pressure and energy density is close to that computed in LQCD. 
Wc have then made estimates of the specific heat, the speed of sound and conformal measure. 

Further, we have extracted the Taylor expansion coefhcients of pressure in the two chemical potentials upto 6-th 
order. All the coefficients approach their respective SB limit above 2Tc. A quantitative comparison with the LQCD 
results show reasonable agreement, though the QNS C2 and the INS C2 on the Lattice are smaller by about 20 %. 
In contrast our earlier estimates [6, 8] of these coefficients without the VdM term showed that C4 and cj differ from 
the LQCD results. Thus the main effect of the VdM term is to impose physical constraints on $ and $ such that 
at large temperatures the coefficients of the same order approach each other. This is clearly visible from the flavour 
off-diagonal coefficients shown in Fig. 7. The remaining difference of the values of the QNS and INS in the model and 
lattice still needs to be addressed. Possible future steps to bring in better agreement could be to include beyond mean 
field effects and/or to include some sort of temperature dependence to the coefficient of the VdM term. However 
the lattice quark masses may be important in bridging the gap. We already found that such data for pressure with 
almost physical quark masses [37] show an increase at any given temperature when compared to data with larger 
quark masses [36] . This would encourage us to believe that extraction of the susceptibilities with similar quark masses 
on the lattice may have a better agreement with our results. Another way to compare results would be to re-estimate 
the parameters of the NJL model directly from the pion mass and decay constants from the lattice. We hope to 
undertake such studies in future. 

In an alternative formulation of the PNJL model including the effect of the VdM term, the coefficients C2, C4, Cg 
and cg have been calculated [11]. Surprisingly, we more or less agree with those results quantitatively. Apart from 
the fact that this may be possible due to various adjustable parameters in both the models, the main reason seems 
to be the small dependence of and $ on the chemical potentials. The basic difference between the two approaches 
is in the use of the VdM potential. The VdM term is required to obtain the mean field solution of $ and $. But as 
we have explained in the formalism that it should not be included in the expression for pressure. On the other hand 
in Ref. [11] apart from obtaining the mean fields the VdM term is included while calculating the value of pressure. 
The difference in the mean field treatment coupled by almost same final results provide hints to the fact that mean 
field treatment has certain shortcomings and is unable to settle issues at hand. It would thus be worthwhile to look 
beyond. 
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APPENDIX A 



We demonstrate the methodology of extracting the average value of a quantity from mean field approximation. 
Suppose wc have a variable ^' with a probability distribution given by cxp(— ^f^) and we have to obtain the average of 
the function in the exponential ^'^ (just like we have to obtain the average of V in Eqn. 2). We define the distribution 
in the domain < 4* < 1. The distribution is shown in Fig. 8(a). The average is obtained as. 



(5»2^ ^ JO ^ 0.2537 (Al) 

where Z = d^e^* = 0.747 is like a partition function. In the given domain the distribution has no maximum and 
thus a mean field solution cannot be obtained. Let us now make a change of variable from \1/ to $ where = <I>^. 
The distribution becomes 2$ exp(— as shown in Fig. 8(b). Here, 2$ is like the Jacobian in the main text. One 
can now easily check that corresponding to Eqn.Al, we need to find the expectation value of given by. 
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(a) (b) 
FIG. 8: (a) $ and (b) pressure for k = 0(ro = 0.27 GeV) and k = 0.5(ro = 0.2555 GeV). The value of Tc is 0.270 GeV. 



($4) = ^° ^^^""^ = 0.2537 (A2) 

In this case the distribution has a maximum and we can do a saddle point approximation. We thus minimize ln[2(f>] 
which gives the mean field value ($) = 1/^2. Using this value we find = 0.25 ~ ($4). So gives a good 
approximation to On the other hand if we include the logarithm term we have. ($)^ — ln[2($)] = —0.0966, 

which is widely different. 
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